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ABSTRACT 

The transiting exoplanet WASP- 18b was discovered in 2008 by the Wide Angle Search for Planets (WASP) 
project. The Spitzer Exoplanet Target of Opportunity Program observed secondary eclipses of WASP- 1 8b using 
Spitzer's Infrared Array Camera (IRAC) in the 3.6 |a.m and 5.8 [xm bands on 2008 December 20, and in the 4.5 
^m and 8.0 \im bands on 2008 December 24. We report eclipse depths of 0.30±0.02%, 0.39 ±0.02%, 0.37 ± 
0.03%, 0.41 ±0.02%, and brightness temperatures of 3100±90, 3310±130, 3080±140 and 3120±110 K in 
order of increasing wavelength. WASP- 18b is one of the hottest planets yet discovered - as hot as an M-class 
star The planet's pressure-temperature profile most likely features a thermal inversion. The observations also 
require WASP- 18b to have near-zero albedo and almost no redistribution of energy from the day-side to the 
night side of the planet. 

Subject headings: planets and satellites: atmospheres — planets and satellites: composition — planets and 
satellites: individual (WASP- 18b) — infrared: planetary systems 



1. INTRODUCTION 

Among the more than 500 extrasolar planets discovered to 
date', the over 100 close-orbiting gas giants that transit their 
host stars have provided the most valuable clues to their phys- 
ical natures. The geometry of the transit gives a direct mea- 
surement of the density of the host star and the surface grav- 
ity of the planet (Seager & Mallen-Ornelas 2003; Southworth 
et al. 2007). These measurements can be combined with an 
estimate of the star's mass and radius to provide estimates 
of the planet's mass, radius and density. The closest-orbiting 
planets attain dayside temperatures high enough to give ob- 
servable secondary eclipses at thermal-infrared wavelengths 
as they pass behind their host stars (e.g., Charbonneau et al. 
2005; Deming et al. 2005; Sing & Lopez-Morales 2009). 

Observations with the Spitzer Space Telescope have re- 
vealed that transiting gas giant planets can be divided into 
two classes based on their infrared spectral energy distribu- 
tions. A subset of the very hottest planets, with dayside tem- 
peratures in excess of 2000 K, display molecular features of 
CO and H2O in emission rather than absorption, indicating the 
presence of a temperature increase with height in the planet's 
photospheric layers. Such temperature inversions have been 
inferred from the flux ratios between the four bandpasses of 
the Infrared Array Camera (IRAC, Fazio et al. 2004) aboard 
Spitzer, centered at 3.6, 4.5, 5.8, and 8.0 |a.m (channels 1-4), 
in the hot planets HD209458b (Burrows et al. 2008; Knut- 
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son et al. 2008), XO-lb (Machalek et al. 2008), TrES-2b 
(O'Donovan et al. 2010), TrES-4b (Knutson et al. 2009a), 
XO-2b (Machalek et al. 2009) and WASP- lb (Wheatley et al. 
2011), among others. This phenomenon has been attributed 
to the presence of strongly-absorbing species such as TiO and 
VO remaining in the gas phase in the upper atmosphere, lead- 
ing to the formation of a stratospheric temperature inversion 
(Burrows et al. 2008; Fortney et al. 2008). This led Fortney 
et al. to propose a scheme in which planets are assigned to 
class pM (with hot stratospheres) or pL (without) according 
to temperature, by analogy with the stellar M and L spectral 
classes. Knutson et al. (2010) suggest instead that some plan- 
ets lack inversions because the TiO and VO are destroyed by 
UV radiation, while Spiegel et al. (2009) claimed that TiO/VO 
may not be sufficiently abundant in the upper atmospheres to 
produce the required inversions. Zahnle et al. (2009b) pro- 
posed additional absorbers, like sulphur species, which could 
potentially form thermal inversions. 

WASP- 18b is a prime candidate for secondary-eclipse ob- 
servations to test for the presence of a hot stratosphere (Hellier 
et al. 2009). It orbits an F6-type star with a period of just 0.94 
days and is expected to attain a blackbody equilibrium tem- 
perature approaching 2400 K, assuming zero albedo and effi- 
cient transport of heat from the dayside to the nightside of the 
planet. With less efficient transport, an even higher dayside 
temperature is expected. The planet is unusual because of its 
very high mass, Mp = 10.43 ± 0.30 Jupiter masses, Mjup, and 
modest radius, Rj,= 1 . 165 ± 0.055 Jupiter radii, Rj^p (South- 
worth et al. 2009), which give it a surface gravity an order 
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TABLE 1 
Run Parameters 



A 
Obs. Date (i^m) 


Pre-Observation 
Start" Frame Time (s) 


Frames 


Main Observation 
Start" Frame Time (s) 


Frames 


Post-Observation 
Start" Frame Time (s) 


Frames 


CH13 2008-12-20 3.6 5.8 
CH24 2008-12-24 4.5 8.0 


820.57796 2, 12 
824.35796 12 


260 
185 


820.61869 2, 12 
824.38431 2, 12 


1148 
1148 


820.79721 2, 12 
824.56278 2, 12 


10 
10 



" BJD - 2,454,000 ephemeris time 

of magnitude greater than that of any known transiting planet 
likely to belong to the pM class. The pressure scale height in 
its photospheric layers should therefore be an order of magni- 
tude smaller than for planets of similar temperature, provid- 
ing an important new dimension for tests of planetary atmo- 
spheric models in the presence of strong external irradiation. 
The small but precisely known orbital eccentricity e=0.0085 
±0.0008 (Triaud et al. 2010) may impart slightly faster-than- 
synchronous rotation, which could help to redistribute heat 
from the dayside to the nightside of the planet. 

WASP- 18b was therefore selected as a candidate for obser- 
vation as part of our Spitzer Exoplanet Target of Opportunity 
Program shortly after its discovery was confirmed. In Sec- 
tions 2 and 3 of this paper we describe the observations and 
the analysis of the data. In Section 4 we discuss the con- 
straints imposed on the thermal structure of the planetary at- 
mosphere by the IRAC fluxes, in Section 5 we compare the 
eccentricity and orientation of the orbit from the timing and 
duration of the secondary eclipse with the values derived from 
radial-velocity observations, and in Section 6 we present our 
conclusions. 

2. OBSERVATIONS 

Spitzer's IRAC instrument observed (program 50517) two 
secondary eclipses (see Table 1). After each observation, an 
offset 6-minute full-array sequence confirmed the lack of per- 
sistent bad pixels near the stellar position. There are two in- 
dependent analyses of these data, one presented here and one 
byMachaleketal. (2011). 

IRAC exhibits some sources of systematic noise that must 
be taken into account when planning and analyzing observa- 
tions. A positional sensitivity exists in the 3.6 and 4.5 \xm 
channels (Charbonneau et al. 2005; Knutson et al. 2009b; 
Machalek et al. 2011; Stevenson et al. 2010), and a time- 
varying sensitivity ("ramp") exists in the 5.8 and 8.0 [a.m chan- 
nels (Charbonneau et al. 2005; Harrington et al. 2007). 

The time-varying sensitivity in the 8.0 [a.m channel mani- 
fests as an apparent increase in flux with time. The rate of 
increase depends on the number of photons received by each 
pixel, and is believed to be caused by charge trapping. The 
effect is successfully reduced by staring at a bright diffuse 
source, in this instance an HII ionized region, prior to the 
main observation (a "preflash"). The large number of pho- 
tons quickly saturates the detector, resulting in a smaller rate 
of increase in the main observation than is seen without a pre- 
flash (Knutson et al. 2009b). We observed a 30-minute pre- 
flash prior to the December 24 event (see Figure 1), which ex- 
hibits a decreasing ramp, unlike previous observations (Knut- 
son et al. 2009b; Campo et al. 2011). This is attributable to 
the previous IRAC observation of the bright extended source 
IC1396a. 

The 5.8 |xm channel exhibits a decreasing ramp, so Spitzer 
stared at the target for 62 minutes prior to the observation to 
allow time for the detector to stabilize. In order to minimize 
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Fig. 1 . — Preflash light curve. These are 8 |j.m data, analyzed with aperture 
photometry. The preflash source is bright compared to WASP- 18, which al- 
lows the array sensitivity to stabilize before the science observations. Without 
a preflash, similar observations generally show a steeper and longer ramp in 
the eclipse observations. This may be the first descending preflash, attributed 
to an even brighter source in the immediately preceeding program. 



the positional sensitivity at 3.6 and 4.5 [a.m, each observation 
used fixed pointing. 

3. DATA ANALYSIS 

All data are Spitzer Basic Calibrated Data (BCD) frames us- 
ing version S18.7.0 of Spitzer's preprocessing pipeline. This 
pipeline removes all well-understood instrumental signatures 
and produces a flux-calibrated image (Fazio et al. 2004). We 
first account for light travel time in the solar system by con- 
verting to Barycentric Julian Date (BJD), and then mask all 
pixels in the Spitze r-supplied permanently bad pixel mask. 
We find the remaining bad pixels by grouping sets of 64 
frames and doing a two-iteration outlier rejection at each pixel 
location. Within each array position in each set, this routine 
calculates the standard deviation from the median, masks any 
pixels with greater than Aa deviation, and repeats this proce- 
dure once. Masked pixels do not participate in the analysis. 

WASP- 18 is very bright relative to the background. A 2D 
Gaussian fit to data within 4 pixels of the stellar brightness 
peak determines the stellar center in each frame. The light 
curve comes from 5 x -interpolated aperture photometry (Har- 
rington et al. 2007), excluding frames with masked pixels in 
the photometry aperture and not using masked pixels in sky 
level averages. Table 2 presents photometry parameters. We 
vary the aperture radius between 2.0 and 5.0 pixels in 0.25- 
pixel increments, choosing the one with the best light-curve 
fit as described below. 

We model the intrapixel variation affecting the 3.6 and 4.5 
l^m channels with a second-order, two-dimensional polyno- 
mial (Knutson et al. 2008; Stevenson et al. 2010; Campo et al. 
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TABLE 2 

Details of the Analysis and Results From Light Curve Fit 



Parameter 


3.6 i^m 


4.5 |.tm 


5.8 |.tm 


8.0 i^m 


Array Position (x, pix) 


30.24 


23.20 


24.41 


24.17 


Array Position (y, pix) 


23.89 


24.68 


24.05 


22.32 


Position Consistency" (Sjc, pix) 


0.018 


0.019 


0.028 


0.022 


Position Consistency" (Sy, pix) 


0.018 


0.038 


0.015 


0.021 


Aperture Size (pix) 


3.50 


3.00 


3.75 


3.50 


Sky Annulus Inner Radius (pix) 


7 


7 


7 


7 


Sky Annulus Outer Radius (pix) 


12 


12 


12 


12 


System Flux (f,, \i3y) 


168080 ±140 


104300 ± 600 


69690 ± 20 


37450 ±10 


Eclipse Deptli (Fp/F„ %) 


0.30 ±0.02 


0.39 ± 0.02 


0.37 ±0.03 


0.41 ±0.02 


Brightness Temperature (K) 


3100 ±90 


3310±130 


3080 ±140 


3120±110 


Eclipse Mid-time (?mid, phase)'' 


0.4995 ±0.0007 


0.4985 ± 0.0006 


0.4995 ± 0.0007 


0.4985 ± 0.0006 


Eclipse Mid-time (r^a, BJD - 2,454,000)'' 


820.7160 ±0.0006 


824.4809 ± 0.0005 


820.7160 ±0.0006 


824.4809 ± 0.0005 


Eclipse Duration (r4_i , sec)*" 


8010 ±60 


8010 ±60 


8010 ±60 


8010 ±60 


Ingress (i2-i) and Egress ((4-3) Times (sec)*" 


857 


857 


857 


857 


Ramp Name 


Linear 


Linear 


Falling Exponential 


Rising Exponential 


Ramp, Linear Term (ri ) 


0.005 ±0.001 


-0.006 ±0.003 








Ramp, Curvature (r2) 








14 ±1 


17±1 


Ramp, Time Offset (rs) 








-0.035689 


0.082618 


Intrapixel, Quadratic Term in x (7)2) 





0.47 ±0.13 








Intrapixel, Cross Term (ps) 





-0.12 ±0.01 








Intrapixel, Linear Term in y (p4) 


0.067 ±0.004 











Intrapixel, Linear Term in x (pj) 


-0.086 ±0.003 


-0.33 ± 0.06 








Total frames 


1148 


1148 


1148 


1148 


Good frames 


1142 


987 


996 


1031 


Rejected frames (%) 





14 


13 


10 


Standard Deviation of Normalized Residuals 


0.002428 


0.003485 


0.003738 


0.002928 


Uncertainty scaling factor 


0.31674 


0.47164 


0.57274 


0.47792 



" RMS frame-to-frame position difference 

^ Duration and ingress/egress time are each a single parameter shared among all four wavelengths. 

parameter for each pair of channels observed together. 



Eclipse mid-times are a single 



2011), 



Vip(x,y) = piy^ + P2X^ + P3xy+ p4y+ P5X+ I, 



(1) 



where x and y are the centroid coordinates relative to the pixel 
center nearest the median position, and pi - ps are free pa- 
rameters. The systematics had little to no dependence on the 
quadratic y term, so /?i is fixed to zero for all models. The 
ramps for the 3.6 and 4.5 [a.m channels use a linear model. 



/?(f) = ri(f-0.5)+l, 



(2) 



where t is orbital phase and r\ is a free parameter. We model 
the 5.8 [a.m channel with a falling exponential 



/?(f)=l+exp(-r2[f-r3]), 



(3) 



and the 8.0 [a.m channel with a rising exponential (Harrington 
et al. 2007) 



/?(f)=l-exp(-r2[f-r3]), 



(4) 



where r2 and r^, are free parameters. The rj, term is fixed to its 
best-fit value. The eclipse, E{t), is a Mandel & Agol (2002) 
model which includes the time of secondary eclipse, t\ to 14 
duration (P' to 4* contact), ingress/egress time, and eclipse 
depth. 
The single-channel light curve model is 



Fix,y,t) = F,Vipix,y)Rit)Eit), 



(5) 



where F(x,y,t) is the flux measured from interpolated aper- 
ture photometry and F^ is the (constant) system flux outside 
of eclipse, including the planet. We dropped a small num- 
ber of initial frames in each light curve (0, 150, 100, and 80 
frames in order of increasing wavelength) to allow the point- 
ing and instrument to stabilize. The model lines in Figure 



2 show which points are included, as do the electronic light 
curve files. We fit Eq. 5 to the data using a least-squares mini- 
mizer. Because the Spitzer pipeline usually overestimates un- 
certainties, we re-scale the photometric uncertainties to pro- 
duce a reduced x^ of 1 and re-run the fit. This typically con- 
verges in one iteration. For a given photometric set the scaling 
factors are almost identical for all models, so we choose one 
for use with all models. 

Because the underlying physics of the systematics have not 
been characterized sufficiently to find an expression that fits 
well in every case, it is possible that, for some Spitzer data 
sets, different investigators will find different values for key 
parameters when using different systematic models and pho- 
tometry parameters. This occurred, for example, when Knut- 
son et al. (2009b) re-analyzed the data of Harrington et al. 
(2007), finding an eclipse about half as deep. The x^ mini- 
mum that Knutson et al. found was present in the correlation 
plots of Harrington et al., but it was just a local minimum for 
that model. While models for deep eclipses, such as those pre- 
sented here and by Campo et al. (201 1), should generally pro- 
duce compatible results even with different systematic mod- 
els, weak eclipses such as those of Harrington et al. (2007) 
and Stevenson et al. (2010) are more dependent on the details 
of fitting. 

All of the published results of our current pipeline (Steven- 
son et al. 2010; Campo et al. 201 1, and this work) result from 
testing a variety of models for each systematic, as well as as- 
sessing the best photometry aperture and stabilization time for 
each data set. We test linear, quadratic, quartic-in-log-time, 
falling or rising exponential, logarithmic-plus-quadratic, and 
logarithmic-plus-linear ramps, and a variety of polynomial in- 
trapixel models, before choosing the final models. Most of the 
possibilities produce obvious bad fits. For this paper we se- 
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Fig. 2. — Raw (left), binned (center), and systematics-corrected (right) secondary-eclipse light curves of WASP-18b in the four IRAC channels, normalized 
to the mean system flux within the fitted data. Colored lines are the best-fit models; black curves omit their eclipse model elements. A few initial points in all 
channels are not fit. as indicated, to allow the telescope and instrument to stabilize. 



lect the best two models in each channel and fit them for all 
apertures. 

For each channel, photometry using the various aperture 
sizes produces slightly different data sets. We must select first 
the best data set and then the best model, but x^ and related 
fitting criteria only compare different models to a single data 
set; they are inappropriate for deciding between models fit to 
different data sets. For data sets from different apertures, we 
choose the one with the smallest standard deviation of nor- 
malized residuals (SDNR) with respect to the system flux for 
a given model and repeat for each of several models. This 
generally results in a consistent choice of the best aperture 
size among different models. 

Once we have the optimal aperture size, we then compare 
the models. Since adding additional parameters to a model 
will always produce a better fit, we use fitting criteria that 
properly penalize the addition of parameters. As described by 
Campo et al. (2011), we apply both the Akaike Information 
Criterion, 



AIC = x +2/t, 



(6) 



where k is the number of free parameters, and the Bayesian 
Information Criterion, 



mC = x +k\nN, 



(7) 



where A^ is the number of data points (Liddle 2007). A lower 
information criterion value indicates a better model. Figures 
3, 4, 5, and 6 present SNDR and BIC with the two main candi- 
date models and aperture sizes for each wavelength. Our final 
joint model fit, with 28 free parameters, combines the eclipse 
durations for all channels and pairs the simultaneously ob- 
served mid-times. It resulted in an AIC of 4176 and a BIC of 
4302. 

To assess parameter uncertainties and correlations we ex- 
plore phase space with a Metropolis random-walk Markov- 



chain Monte Carlo (MCMC) routine. Each chain began at the 
least-squares minimum. If any step in the chain ever beat the 
minimum, it would indicate an even deeper minimum at the 
bottom of the basin of attraction just entered, so the routine 
would discard the MCMC data, re-run the minimizer, and re- 
start the Markov chain. The routine runs a "burn-in" of at 
least 10^ iterations to forget the starting conditions, and then 
runs four million iterations. 

We also consider the level of correlation in the residuals. 
For this, we plot root-mean-squared (RMS) model residuals 
vs. bin size (Pont et al. 2006; Winn et al. 2008) and compare 
to the theoretical 1/v^RMS scaling. Figure 7 demonstrates 
the lack of significant photometric noise correlation in our fi- 
nal models. In the case of Channel 3, we found a high degree 
of correlation between some of the model parameters in the 
posterior distribution, and prefer a less-correlated model with 
insignificantly poorer BIC and similar SDNR at 3.75 rather 
than 4.0 pixel aperture size. Differences in interesting param- 
eter values for the near-optimal alternative are < Icr. 

Finally, the marginal posterior distributions (i.e., the pa- 
rameter histograms) and plots of their pairwise correlations 
help in assessing whether the phase space minimum is global 
and in determining parameter uncertainties. We present these 
plots for the astrophysical parameters in Figures 8, 9, 10, and 
1 1 . Table 2 gives the values and uncertainties of all parame- 
ters. 

The data files containing the light curves, best-fit mod- 
els, centering data, etc., are included as electronic supple- 
ments to this article. Multiple teams analyze the same Spitzer 
exoplanet data, sometimes obtaining divergent results (e.g., 
Beaulieu et al. 201 1 and Stevenson et al. 2010; Knutson et al. 
2009b and Harrington et al. 2007). To facilitate compari- 
son of these efforts, we encourage all investigators to make 
similar disclosure in future reports of exoplanetary transits 
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Fig. 3. — Channel 1. SDNR and ABIC vs. aperture size. A lower value 
indicates a better model fit. 
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Fig. 6. — Channel 4. SDNR and ABIC vs. aperture size. A lower value 
indicates a better model fit. 



0.003675 






' 








♦ — • 


No Ramp SDNR 




♦ ♦ No Ramp ABIC 


0.003650 


- • — • 


Linear Ramp SDNR 






▼ ▼ Linear Ramp ABIC 


0.003625 








X A^ ^~~~^ 


0.003600 








xy ^ 


a: 0.003575 

Q 

'^ 0.003550 


♦ 








0.003525 


-\\ 


, 


^x 




0.003500 


\ 


v-:^ 


♦ 


♦ 
♦ ♦ ♦ T 


0.003475 










0.003450 











y 

< 



2.50 2.75 3.00 3.25 3.50 3.75 4.00 4.25 4.50 
Photometry Aperture Size [Pixels] 

Fig. 4. — Channel 2. SDNR and ABIC vs. aperture size. A lower value 
indicates a better model fit. 
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Fig. 5. — Channel 3. SDNR and ABIC vs. aperture size. A lower value 
indicates a better model fit. 



and eclipses. Because of differing photometry methods and 
the vagaries of estimating error, the standard deviation of the 
residuals, normalized to the out-of-eclipse flux, should be the 
figure of merit for comparing analyses of the same data by 
different pipelines. 







Bin Size 



Bin Size 



Fig. 7. — Root-mean-squared (RMS) residual flux i'.v. bin size in each chan- 
nel. This plot tests for correlated noise. The straight line is the prediction for 
Gaussian white noise. Since the data do not deviate far from the line, the 
effect of correlated noise is minimal. 



4. ATMOSPHERE 

We model WASP- 18b with the exoplanet atmospheric mod- 
eling and retrieval technique developed by Madhusudhan & 
Seager (2009). This is a ID, line-by-line, radiative-transfer 
model with constraints of hydrostatic equilibrium and global 
energy balance. The model has six temperature structure pa- 
rameters and four molecular abundances, expressed as devia- 
tions from thermochemical equilibrium and solar abundances. 
Recognizing the excess of model parameters over data, our 
goal is to rule out unreasonable areas of phase space rather 
than to determine a unique composition and thermal profile. 
An MCMC routine runs a wide range of inversion and non- 
inversion models, integrates the resulting spectra against the 
Spitzer bandpasses, and calculates x^ against the four data 
points. Integrals over the MCMC posterior distribution pro- 
duce robust statistical statements about the underconstrained 
model. 

Figure 12 shows the observed planet-star flux ratios and 
two model spectra. We find that the observations can be ex- 
plained by models with and without thermal inversions. We 
note that the observations are also consistent with a black- 
body planetary spectrum with T = 3200 K, although a black- 
body spectrum is Hkely unrealistic. An atmosphere can have 



Nymeyer et al. 



<= 0.102 

B 0.101 

^ 0.100 

a 0.099 

S 0.098 

= 0.097 

■^ 0.096 

0.0036 

JO 0.0034 

B^ 0.0032 

iij X 0.0030 

^i^ 0.0028 

u 0.0026 

0.0024 

168600 

e ;^. 168400 

^ 3. 168200 

i^x 168000 

^ iZ 167800 

u 167600 

^ 167400 

;S 0.5005 

£ 0.5000 

°- 0.4995 

S 0.4990 

.9- 0.4985 

,y, 0.4980 

■; 0.4975 

S 0.4970 

I 0.4965 

0.0050 

|_| 0.0045 

iij X 0.0040 

Ji^ 0.0035 

0.0030^ 




CH13 Eclipse Phase Eclipse Durati 



CHI Eclipse CHI System CH24 Eclipse Phase 

Flux Ratio Flux (^/y) 




Fig. 8. — Parameter correlations. To decorrelate the Markov chains and 
unclutter the plot, one point appears for evei'y 1000th MCMC step. Each 
panel contains all the points. 



Fig. 10. — Parameter correlations, continued. Same as Figure 8. 
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Fig. 9. — Parameter correlations, continued. Same as Figure 8. 



a blackbody spectrum either if it is isothermal over the upper 
several optical depths or if there is no opacity source. Nei- 
ther condition is physically favorable. Several spectrally ac- 
tive molecules should be abundant in hot- Jupiter atmospheres 
and there is collision-induced opacity (Freedman et al. 2008). 
In addition to being coupled with the opacities, the tempera- 
ture structure is also critically influenced by atmospheric dy- 
namics (Showman et al. 2009), all of which can cause a non- 
isothermal profile. 

At 3200 K, the temperature of early M-class stellar photo- 
spheres, CO and H2O are the dominant spectroscopically ac- 
tive molecules in the IR. Other molecules like CH4 and CO2 
are neghgible, under the assumption of thermochemical equi- 
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Fig. 1 1. — Pai'ameter histograms. To decorrelate the Markov chains, the 
histograms come from every 100th MCMC step. 



librium with solar abundances. CO has a strong absorption 
feature in the 4.5 ^.m channel. H2O contributes the dominant 
opacity in the 5.8 ^m channel, and contributes significantly 
in the remaining channels. Thus, for temperature decreasing 
monotonically with altitude, i.e., in the absence of a thermal 
inversion, the spectra should exhibit noticeable absorption in 
the 4.5 and 5.8 |^m channels and less in the 3.6 |J.m and 8 
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l^m channels (Madhusudhan & Seager 2010, 201 1). The ob- 
served planet-star flux contrast in the 4.5 \im channel should 
then be lower than that in the 3.6 iJ.m channel (Charbonneau 
et al. 2008; Stevenson et al. 2010); the difference depends on 
the temperature gradient and the composition. 

Our observations of WASP- 18b show excess flux at 4.5 |xm, 
compared to the 3.6 ^.m channel. This could be due to a ther- 
mal inversion, but the observational uncertainties also allow 
just a gentle temperature gradient and no inversion, such that 
the absorption features are not too deep, and a different chem- 
ical composition (Madhusudhan & Seager 2010). Two mod- 
els, with and without an inversion, appear in Figure 12. Both 
explain the data fairly well, the inversion model at the Icr 
level and the non-inversion model within 1.5cr. The molec- 
ular abundances of the models are only marginally different 
from those of thermoc hemic al equilibrium with solar abun- 
dances (TEq). The inversion model has 10 times more CO, 
and the non-inversion model has 10 times less H2O and CO, 
as compared to TEq. The mixing ratio of CO2 is 10'^ for the 
inversion model and 10"^ for the non-inversion model. Despite 
the weak constraints on the temperature structure, the obser- 
vations do place a strict constraint on the day-night energy 
redistribution in WASP- 18b: Both models require a low Bond 
albedo (A) and inefficient day-night redistribution (< 0.1 for 
A=0) in WASP- 18b. Figure 13 shows the contribution func- 
tions of the two models in the four IRAC channels, along with 
the thermal profiles. 

The presence of a thermal inversion in the dayside atmo- 
sphere of WASP- 18b is expected based on theoretical grounds 
(Hubeny et al. 2003; Fortney et al. 2008). At the high temper- 
atures of this planet, TiO and VO can exist in gas phase over 
the entire atmosphere, thus contributing to the strong visible 
opacities required to form stratospheres. However, questions 
of whether the concentrations of TiOA'^O alone are adequate 
to cause the required thermal inversion, and whether other 
sources of visible/UV opacity are possible at these temper- 
atures, merit future theoretical investigation (Spiegel et al. 
2009; Zahnle et al. 2009a). Furthermore, Knutson et al. 
(2010) find that the host star WASP- 18 has a low activity level, 
indicating that inversion-causing compounds are not likely to 
be destroyed by stellar UV radiation, thereby also favoring the 
presence of a thermal inversion. Thus, WASP- 18b is an apt 
candidate for follow-up observations searching for thermal in- 
versions. Stronger constraints on the temperature structure of 
WASP- 1 8b are possible in the near future if ground-based ob- 
servations of thermal emission become available (Madhusud- 
han et al. 201 1). Also, the low day-night energy redistribution 
required by the present observations can be verified by poten- 
tial future observations of thermal phase curves (e.g., (Knut- 
son et al. 2008)) of WASP- 18b. 

5. ORBITAL DYNAMICS 

Our secondary eclipse times further constrain the planet's 
already-precise orbital parameters. Triaud et al. (2010) de- 
tect an eccentricity for WASP-18b of e = 0.0085 ± 0.0008, 
the lowest fully-determined value for any transiting planet 
measured with such precision. A joint photometric fit to all 
four Spitzer observations yields a midpoint phase of 0.4990 ± 
0.0004 for the Hellier et al. (2009) ephemeris, and a duration 
of Ds = 0.0927 ±0.0007 days, which is longer than the transit 
duration of HeUier et al. (2009) by 4.8 a. By itself, after a 
20-second light-time correction, the eclipse midpoint tells us 
that ecosw = -0.0016 ± 0.0007, where u! is the longitude of 
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Fig. 12. — Dayside spectrum of WASP- 18b. The black circles with error 
bars show our observations of WASP-18b in the four Spitzer IRAC channels. 
The red curve shows the inversion model spectrum and the green curve shows 
the non-inversion model spectrum discussed in the text. The red and green 
circles are the respective spectra integrated over the Spitzer bandpasses (in- 
dicated with arbitrary scale at the bottom). The black dashed line shows a 
blackbody at 3150 K, and the blue dashed lines show blackbody spectra cor- 
responding to the minimum and maximum temperatures in the atmosphere 
(see Figure 13). 
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Fig. 13. — Pressure-temperature profile (left) and contribution functions 
(middle and right). The middle and right panels show the normalized contri- 
bution functions for the non-inversion and inversion models, respectively, in 
the indicated Spitzer filters, with wavelengths in ^.m. The left panel overlays 
the profiles for both models. 



periastron. We combine the eclipse phase and duration with 
known transit parameters from Hellier et al. (2009) to deter- 
mine that e sin w = 0.0198 ±0.0072. To do this, we use Eq. 18 
from Ragozzine & Wolf (2009) as derived by Kopal (1959): 



esmu} ■■ 



2 2- 

a -COS I 



\Ds + DpJ \a^-2cos^i 



(8) 



Dp is the transit duration and Ds is the secondary eclipse dura- 

R^ 1+esincj 
a \-e- ' 

are the planetary 



tion. We define a = ^(1 -I- ^) . ^ and cos; = b^ 
and ^ 



y^ 



where /?p and /?* 



and stellar radii, respectively, a is the orbit's semi-major axis, 
P is the orbital period, and hf is the impact parameter in units 
of the stellar radius. We solved the equation numerically for 
esinoj, and the uncertainties come from sampling Gaussian 
distributions generated from the uncertainties of the input pa- 
rameters. 

We joindy fit an MCMC orbit model (Campo et al. 2011) 
to the BJD time for 2 eclipses, 37 radial velocity (RV) data 
points from Triaud et al. (2010), and 6 transit midpoints ex- 
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traded from the photometry of Southworth et al. (2009), omit- 
ting three RV points subject to the Rossiter-McLaughlin ef- 
fect. This results in a total of 45 data points, 42 of which 
are included. Our model has 6 free parameters. Results are 
given in Table 3, where Tq is the ephemeris time, K is the RV 
amplitude, and 7 is the barycentric velocity. 

We find from the fit above that ecosu; = 
-0.00014±0.00053, consistent with 0, leaving e domi- 
nated by its esinoj component. We note that the value of 
esinoj found photometrically is positive, while the value in 
our joint fit is negative. Eclipse timing does not effectively 
constrain e sinoj. Radial velocity measurements have a known 
tendency to overestimate e when it is low (Laughlin et al. 
2005). A key sign of this is a value of a; ^ ±90°. Photometric 
information, such as the measured durations for our eclipses, 
cannot yet independently confirm a non-zero esinw beyond 
the 3(7 level. Precise determination of e is important because 
the circularization timescale for a tidal damping quality 
factor, Qp ^ 10*" (Mardling 2007), is comparable to the age 
of the system. Whether the orbit is still eccentric determines 
if it is still experiencing tidal dissipation, which drives the 
evolution of the system. We perform a comparison fit with 
e = 0, but its BIC value of 155 is considerably higher than 
BIC = 102 for the eccentric fit. 

To determine if the eccentricity we found could have come 
up by random chance, we also performed an experiment sim- 
ilar to that of Laughlin et al. (2005) in which we generated 
10^ radial velocity datasets for a planet in a circular orbit with 
the same period, mid-transit time, and semi-amplitude as our 
best fit for WASP- 18b. In each dataset, the BJD of each ob- 
servation was kept the same as in the real WASP- 18b dataset. 
We added Gaussian noise corresponding to the instrumental 
error for each observation and a 3 m s ' stellar jitter consis- 
tent with the tables of Wright (2005), in quadrature. We re- 
tained the transit and eclipse timings. We then used a min- 
imizer to find combinations of esinw and ecosu; that min- 
imize the x^ corresponding to each dataset. The mean and 
standard deviation of esinw and ecosw were 0.0000±0.0012 
and -0.0003 ± 0.0002, respectively. The 3cr upper limit on 
|esina;| is 0.0036, well below our best-fit value in the real 
dataset. This is not surprising because of the high signal-to- 
noise ratio and even sampling (Triaud et al. 2010) of the RV 
data i-r^ -- 180). Given the MCMC results, the BIC com- 
parison between the fits to a circular and non-circular orbit, 
the Monte Carlo experiments, and the weak but consistent 
photometric support, and despite the improbable orientation 
of the orbit, the la non-zero eccentricity of WASP-18b's or- 
bit is likely not an overestimate. 

Despite WASP-18b's short period and close proximity to its 
host star, its high density and low eccentricity make it a poor 
candidate for the detection of apsidal precession (Ragozzine 
& Wolf 2009; Campo et al. 2011). The precession should 
manifest itself as an eclipse/transit timing variation with a pe- 



riod of 600 years and an amplitude of — , or about four min- 
utes (Ragozzine & Wolf 2009). Given the orbit's current ori- 
entation, apsidal motion could be detectable as a ^7 ms differ- 
ence between the best-fit transit and eclipse periods (Gimenez 
& Bastero 1995). This signal would likely be overwhelmed by 
the modulation of the period due to tidal infall (Hellier et al. 
2009), which could be measurable within a few decades. The 
Applegate effect makes very small transit timing variations 
unmeasurable (Watson & Marsh 2010). 

TABLE 3 
Joint Orbital Fits 



Parameter 


Value 


esintj 


0.0091 ±0.0012 


e cos LC 


-0.00014 ± 0.00053 


e 


0.0091 ±0.0012 


w(°) 


-91 ±3 


P (days) 


0.9414518 ±0.0000004 


To (MJD)" 


1084.79363 ± 0.0001 1 


K (ms^l) 


1818±3 


7 (ms'b 


3327 ± 2 


BIC 


102.0 



' MJD = BJD - 2,454,000 (Terrestrial Time) 



6. CONCLUSIONS 

Spitzer observed two secondary eclipses of WASP- 18b us- 
ing all four channels of the IRAC instrument. A blackbody 
model fits the observed brightness temperatures relatively 
well. Slightly better fits to both inversion and non-inversion 
models exist with the inversion model somewhat preferred. 
Because the planet is so much brighter than its predicted equi- 
librium temperature for uniform redistribution, the model re- 
quires near-zero albedo and very low day-night energy redis- 
tribution. The very small scale height makes this atmosphere 
interesting as an extreme example among irradiated planets. 
The addition of secondary eclipse data also improves the or- 
bital parameters, confirming a slight eccentricity. Files con- 
taining the lightcurves, model fits, source centers, and other 
ancillary data appear as electronic attachments to this article. 
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